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Using high-precision quantum Monte Carlo (QMC) simulations within the framework of dynamical mean 
field theory (DMFT), we show that the anisotropic degenerate two-orbital Hubbard model contains two 
consecutive orbital-selective Mott transitions (OSMTs) even in the absence of spin-flip terms and pair- 
hopping processes. In order to reveal the second transition we carefully analyze the low-frequency part of 
the self-energy and the local spectral functions. This paper extends our previous work to lower temperatures. 
We discuss the nature - in particular the order - of both Mott transitions and list various possible extensions. 



Introduction Recently, the experimental discovery of two consecutive orbital-selective metal-insula- 
tor transitions in the effective 3-band system Ca2-a;Sra;Ru04 triggered a cascade of papers also in the 
theoretical literature. The great interest in these "orbital-selective Mott transitions" can easily be under- 
stood. After initial LDA and DMFT band structure 01111' and strong-coupling (5) calculations, it be- 
came clear that the occurence of OSMTs in Ca2_j;Srj;Ru04 is a non-perturbative correlation phenomenon, 
which generalizes the famous concept of the classical "Mott transition" |6| to multi-band systems with 
inequivalent bands. The orbital-selective nature of the transitions also raises fundamental questions, such 
as, e.g., whether phase transitions in two subsystems merge or remain distinct, depending on the way these 
two subsystems are being coupled. 

The microscopic model used in the literature to describe the OSMT phenomenon theoretically is the 
anisotropic degenerate two-orbital Hubbard model with Q El |9l [Tol CI] [T2| or without fu\ [HI [Tsl [16] 
spin-flip terms and pair-hopping in the Hamiltonian. For the model without these additional terms (the 
" J^-model") it was first believed I13II14I that only one transition exists at low temperatures, i.e., that the 
transitions which occur in the separate orbitals merge upon unison. Since two transitions clearly exist for 
the model with spin-flips and pair-hopping (the " J-model") it was believed 1 8 1 that these additional terms 
cause the two transitions of the subsystems to remain distinct and, hence, make a fundamental difference. 
Here we show, extending previous work in flS], that the J2 -model does, in fact, also contain two distinct 
OSMTs. In this manner we establish the J^-model as a minimal model for the description of the OSMT 
phenomenon. 

The Hamiltonian of the -model. 



describes hopping between nearest-neighbor sites i, j with amplitude t„j for orbital m G {1, 2}, intra- and 
wferorbital Coulomb repulsion parametrized by U and U' , respectively, and Ising-type Hund's exchange 
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Fig. 1 Low-frequency analysis of the self-energy. Main 
panel: weight of singularity at a; = for narrow/wide 
band (filled/open symbols) as extracted from polynomial 
fits to the product ijj'E{uj) [shown in inset for T = 1/32 
and U = 2.8, 2.6, 2.4, 2.2, 2.0, 1.8 (top to bottom)]. 



Fig. 2 Quasiparticle weight Z = m/m* . Inset: Z for 
wide/narrow band (open/filled symbols) across OSMTs 
at Uci « 2.0 and Uc2 « 2.5 for T = 1/40. Main panel: 
detailed study of first transition at Uci reveals hysteresis 
(i.e. a l*^' order transition) only for T < T* ^ 0.02. 



coupling Jz', n^^^ = cj^na^ima ^pi" {T7 i}- As usual, we consider this model with U' = U — 2Jz. 
Specifically, our DMFT (QMC) results for this model, to be presented below, were calculated using semi- 
elliptic densities of states with full bandwidths Wi — 2, W2 = 4, for the "narrow" and "wide" band, 
respectively, and interaction parameters = U/i, U' — U/2. 

It is important to emphasize that high-precision in the QMC simulations is essential for the investigation 
of OSMTs, since the signals of in particular the second transition are rather subtle and can easily be missed. 
In 1 16 1 we showed that our QMC results for, e.g., quasi-particle weights, have relative errors of 0(10"^) 
even near the critical interactions Ud — 2.0 and Uc2 — 2.5 for the first and second transition, respectively. 
This is an improvement of up to two orders of magnitude compared to previous QMC results fTTT^l. The 
high precision in our simulations is (at least in part) due to the reduction of the discretization error achieved 
by complementing QMC data with a high-frequency expansion of the self-energy I17II18I . 

Results In the following, we present QMC results for the -model with band widths and interaction 
parameters as stated before. We briefly discuss one of the primary criteria used in Ref. |[15| for detecting 
the two OSMTs in Fig.n(cf. also inset of Fig.|3} before turning to new results. 

As shown in Fig. ^ a low-frequency analysis of the self-energies associated with both bands clearly 
reveals the orbital-selective character of the first transition at Ud ~ 2.0 (associated with the narrow band) 
and the existence of a second transition (associated with the wide band) at Uc2 ~ 2.5: for U « Ud, 
only the narrow band develops a singularity the weight of which (filled symbols) increases approximately 
linearly with U. In contrast, the corresponding weight for the narrow band remains zero (numerically) up to 
U « Uc2', above this point, again, a linear increase with U is observed. Note that temperature dependencies 
are hardly visible; however, the kink at Uc2 becomes slightly sharper with decreasing T. The extraction of 
the weights shown in the main panel from fits to products ujY.{uj) is illustrated in the inset to Fig.^ 

Metal-insulator transitions may also be identified from discrete estimates of the quasiparticle weight 
Z « [1 — Im T^{iTTT) / (ttT)]^^ (even though this observable is somewhat problematic 1 15 16|): kinks in 
Z, most markedly for the narrow band, reveal a transition at Ud ~ 2.0 in the inset of Fig.|2| a second 
transition for Uc2 ~ 2.5 would become apparent at higher resolution. The main panel of Fig. |2| shows a 
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close-up of Z near \Jc\ for a range of temperatures. Evidently, the observables vary continuously with 
U for T = 1/32 and T = 1/40; only for T = 1/50 = 0.02, we find first ti-aces of stable coexistence, 
signaling a first-order transition. Thus, we estimate T*^ w 0.02 which agrees with recent independent exact 
diagonalization (ED) results ['221, but is much lower than earlier QMC and ED estimates (of T*^ « 0.038 
LMJ and T*i > 0.03 [231, respectively). 

Valuable insight into the nature of phases and transitions can be gained from the spectral functions, 
shown for the narrow and wide band in the left and right hand panel of Fig.|3] respectively. The narrow- 
band spectra show typical Fermi-liquid behavior: the quasiparticle peak becomes narrow before it decays 
and, finally, a gap appears for U > 2.05. In contrast, the wide-band quasiparticle peak remains wide, but 
develops a dip near the Fermi energy before a full gap appears at [/ w 2.6. The density of states at the 
Fermi energy, shown in the inset of Fig.|3] further illustrates the OSMT scenario: At Ud w 2.0 the narrow 
band undergoes a metal-insulator transition while the wide band is hardly affected; the latter becomes 
insulating only above Uc2 ~ 2.5. 

Summary and Discussion We have shown that, contrary to previous statements in the literature, the two- 
band Hubbard model with distinct band widths W2 = 2Wi and interaction parameters U, U' — C//2, = 
/7/4 does in fact describe the occurrence of orbital-selective Mott transitions, as seen experimentally in 
the Ca2-2;Srj;Ru04-system. This "Jz-model" can therefore be considered as a minimal model for the 
theoretical description of OSMTs. We also showed that the method used in this paper, viz. high-precision 
QMC calculations at finite temperatures, is well-suited for the study of OSMTs, provided the correct high- 
frequency behavior of the self-energy is implemented with care and full convergence in the DMFT cycle 
is established. The observables studied in this paper, revealing two consecutive Mott transitions, are the 
spectral functions, the low-energy behavior of the self-energy, and the quasiparticle weights. 

The question concerning the order of the transitions at Uci — 2.0 and Uc.2 — 2.5 is both interesting 
and important. The results presented here suggest that physical quantities, like the quasiparticle weights, 
display a jump at sufficiently low temperatures for U = Ud, so that the first transition is likely of first order 
below T*i K, 0.02. In contrast, these observables are continuous for the J^-model at [/ = Uc2, as expected 
for a second-order transition. At present it would be premature, however, to make definite statements 
concerning the second transition, because high-precision QMC calculations are at present difficult (since 
computationally costly) in the relevant low-temperature regime and also because it is numerically difficult 
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to distinguish a real second-order transition at T > from a narrow cross-over, eradiating possibly from 
SiT — quantum critical point. In spite of these numerical uncertainties, we wish to point out that the 
occurrence of a first-order transition at Ud and a second-order transition at Uc2 are physically plausible. 
A first-order transition for the narrow band at Ud is exactly what one expects on the basis of experience 
with the Mott transition in the single-band Hubbard model |6 19 20 181. The orbital-selective phase 
between Ud and Uc2 is characterized by itinerant electrons in the wide band, interacting with immobile 
(localized) electrons in the narrow band, which is essentially Falicov-Kimball physics. Since the metal- 
insulator transition in the symmetric spinless Falicov-Kimball model is of second order (Hubbard-III-like), 
the same scenario could apply to the -model at Uc2- Our results further show non-Fermi-liquid behavior 
of the wide band for Ud ^ U < Uc2, which is again to be expected on the basis of the FaUcov-Kimball 
analogy. The connection to the Falicov-Kimball model was recently also pointed out by Biermann et al. 
L21J : for a discussion of non-Fermi-liquid physics cf. also 1.22,1 . 

As an outlook, we discuss several extensions of the present work. Clearly, since the transition tem- 
perature at Ud appears to be rather low (T*^ « 0.02), it is important to obtain results at even lower 
temperatures than were considered in this paper. The inclusion of spin flips remains an important problem. 
Since the experimental system Ca2-i^Sr2:Ru04 has three bands, two of which are physically equivalent, it 
may be worthwhile to study the OSMTs in a real 3-band-model. For comparison with the experimental 
situation, additional hybridization between the bands needs to be taken into account; first results 1 10 12J 
suggest that hybridization influences the OSMT dramatically. Since anisotropy is by definition important 
for the OSMT, it is only consistent to take also orbital-dependent interaction strengths into account. Since 
the experimental system is an antiferromagnetic insulator for small doping {x < 0.2), it is important to ex- 
tend the present work to include also magnetic phases. In conclusion, the study of OSMTs will doubtlessly 
remain exciting, yielding unexpected important results in the months to come. 
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